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Abstract. The Kubo fluctuation-dissipation theorem relates the current 
fluctuations of a system in an equilibrium state with the linear AC-conductance. 
This theorem holds also out of equilibrium provided that the system is in a 
stationary state and that the linear conductance is replaced by the (dynamic) 
conductance with respect to the non equilibrium state. We provide a simple 
proof for that statement and then apply it in two cases. 

We first show that in an excess noise measurement at zero temperature, in 
which the impedance matching is maintained while driving a mesoscopic sample 
out of equilibrium, it is the nonsymmetrized noise power spectrum which is 
measured, even if the bare measurement, i.e. without extracting the excess part 
of the noise, obtains the symmetrized noise. 

As a second application we derive a commutation relation for the two com- 
ponents of fermionic or bosonic currents which holds in every stationary state 
and which is a generalization of the one valid only for bosonic currents. As 
is usually the case, such a commutation relation can be used e.g. to derive 
Hcisenberg uncertainty relationships among these current components. 

1 Introduction: Definitions, Kubo formula 

Consider first a system in an equilibrium state and its current operator defined 



e is the charge of each of the particles in the system, its momentum and Xi 
its position. We consider the current in a cube of volume L 3 . For simplicity we 
take L 3 to be a unit cube, L = 1, and write our formulae in one dimension^. 

The current fluctuations are often described by the noise power spectrum 
defined by 



by: 
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where 






(3) 



(...) denotes averaging with respect to a stationary state: 



(A) = Tr Po A (4) 

where p$ is a time independent density matrix: 

p = 0, [H,po] = (5) 

and H is the Hamiltonian. 

Suppose now that the system is driven out of equilibrium by applying an 
AC electric field 

E{t)=Ee iu>t . (6) 

For a weak enough E(t) (this regime is called the linear response regime) the 
resulting current will be of the form 

I(t) = a(iu)Ee luJt (7) 

where <t{lo) is time independent. We define the conductance of the system, 
G d (oj), by 

G d {u) = Rea{u). (8) 

In more general situations, one can perturb the system by various external fields 
and measure other properties beside the electrical current. In such cases Gd{uj), 
(or <r(w)) is called the linear-response coefficient. Here we shall focus on the 
electrical current but our discussion is extendable to the general case. 

Let S eq (cu) denote S(u) of a system in equilibrium. In 1956 Kubo 
derived a fluctuation-dissipation theorem which relates Gd(u>), and S eq (uj). It is 
the following: 

S eq (-uj) - S eq (uo) = 2hioG d {uj). (9) 

Justifiably, Kubo called it a fluctuation-dissipation relation since it relates the 
dissipative properties of the system, Gd{uS), with its equilibrium fluctuations. 
There exists also another relation of this type which was derived by Callen and 
Wclton 5 in 1951, and which is widely known as the fluctuation-dissipation 
theorem. It is: 

\{S eq {-u) + S eq (co)) = 2G d H(^ + e ^ T _ i ) (10) 

where T is the temperature. On one hand, the Callen- Welton relation is valid 
only for a system in equilibrium. On the other hand, in his work Kubo stressed 
that Eq|5| enables a prediction of a nonequilibrium property such as the con- 
ductance by a calculation of an equilibrium one joj (Although he did not rule 
out generalization to nonequilibrium). These are probably the two main rea- 
sons why it is often believed that Eq|U| is not valid for nonequilibrium states. 



However, Eq|5]zs valid in any noncquilibrium state 7,8,9 provided that this 
state is stationary. That is, 

S(-w)-S(w) = 2huG d (w). (11) 

Here S(u>) is given by Eq[31at any stationary noncquilibrium state (i.e. with the 
condition [SJ). Gd(to) is the response with respect to a small perturbation which 
is applied to the system which is already driven out of equilibrium by another, 
not necessarily small, perturbation. Like EqO Eq^Jhc-hds also for interacting 
systems. 

For example, consider a mesoscopic system at zero temperature which is 
driven out of equilibrium by an external DC field. As a result a DC current 
arises. S(u>) will then be the nonsymmetrized shot-noise spectrum related to this 
current. Suppose now that an additional small "tickling" AC field E(t) — Ee tut 
is applied on top of the DC one. As a result also an additional current appears: 

AI(t) = (I(t)) E>0 - (I) E=0 = a(u)Ee iut (12) 

and now 

G d (uS) = Rea(u), (13) 

that is, Gd{u) will then be the linear coefficient relating the new field with the 
new current (it is therefore perhaps more appropriate to call it the differential 
AC-conductance to distinguish it from the one valid when the AC field is applied 
in equilibrium). Eq^J relates the shot-noise spectra and this differential AC- 
conductance. 

In section |21 we give a simple, self contained, derivation of Eq^] In section 
[3] it is shown that noise measurement setup that measures the symmetrized 
power spectrum at zero temperature yields the nonsymmetrized one when used 
in an excess noise measurement (provided that the system-setup impedance 
matching is kept constant while the system is driven out of equilibrium) . This 
is a direct consequence of Eq^] In section 0] Eq^] is used to generalize the 
canonical commutation relations valid for a current in a boson field to the case 
of a fermionic one, provided the commutator is replaced by its expectation value 
in a stationary state. 

2 Derivation of the nonequilibrium Kubo for- 
mula 

Eq^Jwas obtained in ref. [5] by calculating the net absorption from a classical 
EM field and using the relation between this dissipation and the conductance. 
A more mathematical proof is given in Ref. [7]. Here we present a simple and 
systematic derivation which follows closely the original one given by Kubo [3], 
except that, we do not make use of the specific form of the density matrix has 



in equilibrium but only assume it to be time independent as in any stationary 
state. Consider a system described by the Hamiltonian 



n pi 

Ho = E± + V (14) 



2m 

i=l 



where V — V(x\, ..,x n ). To describe the application of a small external alter- 
nating electrical field we rewrite it as usual as 

g = ^ (P,- eJ 4(x„ t ))* 

A — ' 2m 

»=l 

(throughout this section we take c = 1). The scalar potential does not appear 
since we are using the transverse gauge. In this gauge one has 

E(x,t) = -A(x,t). (16) 

We write A(xi,t) = J dxA(x, t)8(xi —x) (so now A(x, t) is no longer an operator 
since Xi is in the <5-function) and keep only first order term in A. H becomes 

" p2 p 

H = H^ + V - dxA ^ wo*) 

i— 1 

where j(x) is given by Eq|21 We now assume that A(x,t) is constant within 
the cube L 3 , and vanishes outside of it. Adding another part of A(x,t), sim- 
ply results in adding the linear response to it, so this assumption is not a re- 
strictive one. It is needed only because we are looking for the conductance 
which is related to I(t) given by Eq^ With the above assumption we have 
J dxA(x,t)j(x) = A{t) f L sj(x) and hence (recall L = 1): 

H = H a - A{t)I. (18) 

Since we are looking for a relation for a single frequency u> we consider the case 
in which A(t) is of the form: 



and thus, by Ea lltil 
By Eq[Hwe have 1 : 



Alt) = Ae wt (19) 
E{t) = -iiuAe luJt . (20) 
AI(t) = -d(uj)iujAe luJt . (21) 



1 I appearing in Eg 1211 is the same as in Eq0 i.e., it is the average of j(x), Eq|5] In the 
presence of the vector potential the proper (gauge-invariant) current is given by ja(x) = 
j(x) — ^p(x)A(t) where p(x) = ~^2 t 8(x — xi) is the density and thus one should replace / by 
Ia — I ~ ~Q^-(t) where Q is the total charge. Since the extra (diamagnetic) term — QA(t) 
is linear in A it may affect the linear response. However, by Eq |2UI ^-QA(t) = i^jQE(t). 
Because of the i in front of the real coefficient, — — Q, this term contributes only to the out-of- 
phase (non-dissipative) part of the current. In other words, defining (in analogy with Eq |12l : 
A/yt(t) = 5- y i(a;)i?e 1 "*, one obtains from all the above ReaA(^) = Recr(u)). Only the real part 
of ir(uj), Gd{u), appears in Ea llll and therefore our use of / instead of Ia is justified. 



and also 

AI(t) = TrSp(t)I (22) 

where 

p(t) -po = S P (t) ~ 0(A) (23) 

is the change in the density matrix due to switching on the perturbation. The 
equation of motion of the density matrix is 

p(t) = ~[H,p(t)}. (24) 

Recalling that the same equation holds for p , that Sp(t) ~ 0(A), and keeping 
only first order terms in A we get 

Sp(t) = - l -[H^5p{t)] + ~A(t)[I,p ]. (25) 

Eq[3B]is solved by substituting Sp(t) — e~i Hnt a(t)ei Hot , solving for a(t), ex- 
pressing the result in terms of Sp(t) and then using Eq0 One obtains 

Sp(t) = -J drA(t - t)[I(-t), po] (26) 

where 

I(t) = e -k Hot Ie-i Hot (27) 

is the Heisenberg current operator of the unperturbed system. Inserting this 
result into Eq^land using TrABC = TrCAB one gets 

M{t)= h] Q d T A (t-T)([I(P),I(-T)]). (28) 

Comparing with Eg 12 II and using Eg fT§1 yields 

poo 

hwa(u) = - dTe~ luJT ([I(0),I(~T)}). (29) 



Finally, taking the real part of the last equation while using v = I and the fact 
that in a stationary state one has (I(0)I(— r)) = (I(t)I(O)), we obtain 



S(-lu) - S(lo) = I dTe luJT ([I( T ),I(0)}) = 2hwG d (cu) (30) 
as in Eciim 



3 First application. Excess noise measurement 



Consider a mesoscopic system at zero temperature coupled to a detection setup 
(also at zero temperature) which is designed in such a way that it will measure 
(as is very often assumed) the symmetrized noise spectrum: 

S m (u) = S syrn (uj) = + S(-lu)). (31) 

S m (u>) stands for the measured spectrum. Such a setup may resemble, for 
example, the one used in ^U]. In JT] it was shown that for a very broad 
class of setups, if one subtracts the noise measured at equilibrium from the 
noncquilibrium one the resulting spectrum will be given by Eq[3| i.e. it will 
be nonsymmetrized 2 . The main assumption used was that the conductance 
remains approximately unchanged while the system is driven out of equilibrium 
so that the latter remains impedance-matched to the detector. This ensures 
that all the extra power emitted by the shot-noise is detected. Claiming that 
such a measurement yields S(lu) may seem to contradict the assumption EgEm 
however we shall now show that there is no inconsistency: Also in the case of 
EqEJthe excess measurement yields S(ui) and not S sym (uj). 
By its definition the measured excess noise is 

S m ,excess{^) = S m (uj) — S m . eq {uj) (32) 

where S mjeq (us) is the noise measured in equilibrium. Assuming Eg 1311 we have 

S m , excess M = ^(S(u) + S(~Uj)) - -(S e q(u) + S eg (-(j)). (33) 

This can be written as 

S m , excess M = S(u) - S eq (uo) + ^-(S(-Ul) - S(ul)) ~ ^(S eq (-Ul) - S eq (uj)). (34) 

Applying Eq^Jwe get 

S m ,excess(u) = S(uj) - S eq (uj) + tuj(G d (uj) - G d {u})). (35) 

Since we assumed that the conductance remains the same in and out of equilib- 
rium, the last term on the right vanishes, and one is left with 

S m ,excess(u) = S(u) ~ S eq (u>). (36) 

Finally, since S(uj) is the emission spectrum |13| and since in equilibrium at zero 
temperature there is no emission, one has 

S eq (uj) = k B T = (37) 

2 This was shown also for setups that include an amplification stage (as in IIUO . which 
is usually the one determining the measured quantity (see e.g., 1181 ). For an analysis of a 
detection without amplification see Refs. 1121 and 1131 . 



and thus 



S m , exc ess (w) = S(u) k B T = (38) 

as asserted in Thus, also in the specific case of Eg 1311 the excess noise 

measurement yields the general results Eqs. [SHI an d EH1 We emphasize that 
contrary to a common view in the literature |14| . Ea l31l is merely a specific case 
and not a general rule. For a concrete example where Eq|^]does not hold see 
Ref.HSJ. 



4 Second application. Commutation relations 
for fermionic current components 

We now apply Ea llll in order to obtain commutation relation for fermionic 
current components. Let us define 

I{uj) = —= / dte lut I(t). (39) 



V2^ J- 

I(t) is hermitian and therefore 

7 t (w)=J(-cj). (40) 
In any stationary state one has: 

(JM> = -= / dte^(I(t)) 

V ^7T J-oo 

= dte wt Vp ie <Bi '(i|/|i)e- iBit = J(w)V2tt(J), (41) 

V27T J-oo : 

where Pi is the probability to be in the eigenstate \i), and 
</(u,)I(u/)) = 

— I / dtdt l e luJt+lu) ' t '^P l e lE ^ t {i\Ie- lH ^- t 'h\i)e- lE ^ . (42) 

Defining 

T+=^(t + t'), T_=t'-t 

1 / 1 

t = T+ + -T-t =T+ - -T_, 



one has 

(JH/(u/)) 



i pOO poo 

— / dT„dT + ^P^\Ie lHT - I\i)e~ lE ^ T - e^+^'^+e 1 ^^'^- 

^ J -oo J -OO > 

/OO 
dr_ (Ie iHT - Ie-' lHT -)S(oj + uj^e'^-^- 
-oo 

pOO 

= 5(u + u') / dr_e-^ T -(/(0)/(T_)). (43) 



Thus, 

(I(u)I(uj')) = S(u + w')S(-w). (44) 

The averaged cosine and sine components of a current, I cs (Cl) and J srl (£1) , 
are Hermitian operators defined by 

I cs (Cl) = I(Cl) + P(Cl), I sn (Cl) = -i{I(Q)-P(Cl)), (45) 

where we defined averaging over a frequency bandwidth by 

X(Cl) = [ dioX{uo) B : [CI- -A,0+ ^A]. (46) 

J b 2 2 

In a current carried by a boson field, I cs (Ci) and I sn {Cl) form a canonical pair 
similar to x and p of an harmonic oscillator (and in that case I(Cl) is analogous 
to the annihilation operator of an harmonic oscillator). That is, 

[I cs (Cl),I sn (Cl)]=if(Cl) (47) 

where f(Cl) is a real c-number which may depend on CI. An example for such a 
case is the current field in an ideal transmission line ^B] , ^7] . EqEI allows one 
to derive uncertainty relations involving the current components which have 
important consequences in the theory of quantum amplification |18 | |19 | .|2(J | . 
However, this equation is generally not valid for a current carried by fermions, 
in which case the above commutator is in general an operator. 

To overcome this problem we shall use Eq^J Inserting Eq|3!5] into [3] inte- 
grating, averaging over the band width B, and making use of Eas 1401 and 1441 one 
gets 

5(0) = (P(Cl)I(Cl)) 
S(-Cl) = (I(Cl)P(Cl)). (48) 

Subtracting these two equations and making use of Ea llll result in 

([I{Cl),P(Cl)}) = 2Ml(5 d (n)A, (49) 



where for simplicity we assume that A is small enough so that Gd(Q) remains 
constant in it. From Eg 1451 one sees 



[i«(n) 



,7 sn (fi)] = 2i[I(n),jt(fi)]. 



(50) 



Combining the last two equations we finally get 



sn 



{ft)}) = 4imG d (Q)A. 



(51) 



Thus, we have transformed Ea llll into the form of commutation relations which 
is valid for any current in a stationary state, whether it is carried by fermions 
or bosons. The usefulness of Eg 15 II stems from the fact that in many cases, the 
conductance Gd(Sl) is the same in a large set of stationary states (as e.g. was 
the case in Sec|3J and therefore, within such a set, I cs {^) and I sn (D,) posses 
properties of an ordinary pair of canonical variables. 

5 Summary and conclusions 

Kubo's fluctuation dissipation theorem holds also outside of equilibrium, as 
long as the system is in a stationary state. As a consequence, excess noise 
measurement of, for example, the symmetrized noise spectrum yields the non- 
symmetrized one. Another consequence is that although the commutator of the 
two components of fermionic currents is not in general a purely imaginary con- 
stant c- number (unlike for their bosonic counterparts), the projection of their 
commutator onto all stationary states having the same conductance, is such a 
c-number. This can be shown to result in Heisenberg constraints on the perfor- 
mance of quantum transistors |21j . 
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